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Abstract 

We study the production, spectrum and detectability of gravitational waves in models 
of the early Universe where first order phase transitions occur during inflation. 

We consider all relevant sources. The self-consistency of the scenario strongly affects 
the features of the waves. The spectrum appears to be mainly sourced by collisions of 
bubble of the new phases, while plasma dynamics (turbulence) and the primordial gauge 
fields connected to the physics of the transitions are generally subdominant. 

The amplitude and frequency dependence of the spectrum for modes that exit the 
horizon during inflation are different from those of the waves produced by quantum vac- 
uum oscillations of the metric or by first order phase transitions not occurring during 
inflation. 

A not too large number of slow (but still successful) phase transitions can leave de- 
tectable marks in the CMBR, but the signal weakens rapidly for faster transitions. When 
the number of phase transitions is instead large, the primordial gravitational waves can 
be observed both in the CMBR or with LISA (but in this case only marginally, for the 
slowest transitions) and especially with DECIGO. 

We also discuss the nucleosynthesis bound and the constraints it places on the param- 
eters of the models. 

PACS number: 98.80.Cq Keywords: inflation, gravity waves, first order phase transitions. 
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A Appendices 

1 Introduction 

Gravitational waves carry valuable information about the physics that produced them, as they 
decouple quite soon from their surrounding. In particular, waves generated during inflation 
could open open up important opportunities to study the early Universe. 

Beside the ever present generation in vacuum via quantum fluctuations, gravitational waves 
can also be sourced by the anisotropic stress tensor of fields and fluids. More specifically, 
important sources are expected to be present when first order phase transitions occur. In the 
literature, there has been great interest in the gravitational waves generated by this kind of 
transitions (see [IHZ])- The analysis, however, has been mostly concerned with transitions such 
as the electroweak oneS [I] or during preheating [2]. 

In this work we instead investigate the production, features and detectability of gravitational 
waves from first order phase transitions during inflation. 

Inflationary models exhibiting this kind of transitions have existed since the early times of 
inflationary theory: for example Guth's Old Inflation was indeed driven by a first order phase 
transition. However, the motivation for this analysis is even stronger today, because of the 
appearance of many metastable vacua in effective theories of gravity and high energy physics, 
where tunnelings and transitions among vacua are expected to occur, possibly during inflation. 
For example, it has been shown that long series of connected minima can exist in the string 
theory landscape [8]. 

Furthermore, recent investigations have fully analyzed [9] a cosmological model alternative 
to slow- roll/chaotic inflation, where the inflationary dynamics is actually driven by several first 
order phase transitions: chain inflation. 

Although expected to be important, the signatures in gravitational waves from these sce- 
narios and models have not been studied: this work intends to fill this gap3. 

1 Which however does not appear to be a first order transition for the present lower bound on the Higgs mass. 

2 A partial analysis of waves emission in the specific setup of Fourth Order Gravity (FOG) was done in [TO] . 
Those results and approach arc different from ours. In an analysis of gravitational waves in a specific 
realization of chain inflation was attempted. Our results are very different from those, as in the analysis 
appears not to be consistent. The work is not published on a scientific journal. 
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This paper is organized as follows: in section [2] we present the setup and the approach we 
follow. The detailed investigation is then divided in sections [31 El El We start in section El 
by studying our setup in terms of the relevant physical parameters and the bounds on them 
coming from the request of self-consistency of the scenario and its description. These bounds 
will affect also the emission and characteristics of the gravitational waves. 

In section H] we analyze the production of the gravitational waves and their features. We 
study various possible sources related to first order phase transitions during inflation. More 
specifically, in section WM we compute the spectrum of waves emitted by the collisions of bubbles, 
which appear to be the strongest source. In section H~5l we study the remaining ones. 

Finally, in section [5] we discuss the detectability of the waves in the CMBR and at interfer- 
ometers. 

We conclude in section [61 The appendices contain useful accessory material. 

2 Setup and approach 

Consider a period of inflation in the early Universe where some first order phase transitions 
occur. The theory describing this scenario could be very complicated, with a potential ex- 
hibiting many metastable minima at different energies, a large number of fields and a complex 
dynamics, with rolling, tunneling and jumping phases as the fields pass through the minima. 

The inflationary dynamics could occur in various ways depending on the behavior of the 
fields: for example via the mechanism of chain inflation [9], or when one scalar field undergoes 
slow-rolling while others tunnel through the minima. 

The first order phase transitions take place via nucleation of bubbles of the new phases 
within the old ones. With the expansion of the bubbles and their collisions, the latent heat 
of the transitions^ is released and converted in a radiation-dominated fluid. Many sources of 
gravitational waves become active due to these dynamics. 

We will simplify the description of this setup. In fact, knowledge of the details of the field 
theory is not necessary for the kind of analysis we are going to perform. Inflation and the phase 
transitions can be described by a series of physical parameters (for example, the time-scales of 
the transitions, the Hubble parameter, the nucleation size of bubbles and a few others to be 
introduced in section [3]). We will study the production and features of the gravitational waves 
using these parameters, without resorting to the complicated field theory description. 

The analysis is nevertheless quite complex: instead of computing the physical parameters 
from first principle via the field theory, we will constrain them on the basis of the requirements 
of consistency of the scenario. Indeed, the phase transitions must not backreact too strongly 
on the background, if we want inflation not to be stopped and to be efficient. At the same time 
they must be successful, reaching percolation and large scale thermalization, in spite of the fast 
expansion of the Universe. 

All of this also affects the physical parameters of the sources of gravitational waves, con- 
straining the amplitude and features of the latter. The constraints are so binding, that we 
will be able, using our analysis, to fully determine if a series of first order phase transitions 
compatible with inflation can produce a detectable spectrum of gravitational waves. 

3 In the following, with the term transition we will always intend a first order one, even if we do not explicitly 
write it to avoid repetitions. 
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An advantage of the approach we use is that the results are of a general nature and can 
be adapted to many different field theory models. Indeed, once a particular model is chosen, 
one can compute from first principles the physical parameters, which will now depend on the 
couplings of the fields and the dynamics of the model, and specialize our results on gravitational 
wave emission to the case of interest, being able to test it. We will do this, in the end, using 
the example of chain inflation. 

3 Analysis of the setup: parameters and bounds 

We need to consider a limited number of physical parameters to analyze the emission and 
features of the gravitational waves. We also need to investigate which of their values allow 
successful transitions together with efficient inflation, homogeneity and isotropy. 

In the following we will use units for which c = h = 1. We also use t for the cosmic time, rj 
for the conformal and define Mp lanck = (87rG) _1 . The '(') indicates derivative by t(rj). 

3.1 Evolution of the background 

We need only two parameters: 

• The Hubble scale H — -, here written in cosmic time (a is the scale factor). 

• The parameter e = — j| indicating the time evolution of the Hubble scale. We consider 
quasi-de Sitter scenarios, where e < 1 during inflation. 

These parameters come from using a Friedman- Robertson- Walker metric for the background 
(see appendix lA.lj) . The latter is appropriate if there are homogeneity and isotropy at least 
above certain scales, that is if the bubble of the phase transitions do reach percolation and large 
scale thermalization, and if the radiation-fluid generated by the collisions also thermalizes. We 
are soon going to discuss these points. 

3.2 Phase transitions 

We take into account the possibility that more than one phase transition occur, by indicating 
each of them with a progressive integer 1 < n < N; N is kept generic. 

We list here the relevant physical parameters, postponing the discussion of their bounds to 
section 13.2.11 

• the decay rate T n per unit time between phases n and n — 1 (also called nucleation, or 
tunneling rate). 

T n is related to the decay rate per unit time and volume T n , which is the quantity usually 
obtained in a field theory model via tunneling action (or free energy if the temperature 
is important) [T2]: T n = f dVptys^n, where Vphys is the physical volume. 

• The time-scale of the phase transition n — > n — 1. 

is the lapse of cosmic time during which most of the bubble nucleate, collide and 
thermalize. In appendix \K.2\ we show its relation with the decay rate and the tunneling 
action. The time-scale in conformal time is indicated with and defined in IA.2I 
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• the energy density Ae n released by the transition among phases n,n — 1. 

In each transition some energy is liberated. The energy density at disposal is Ae n = 
£n ~ £n-i ,where e m indicates the energy density in the phase m. It is carried by the 
bubble walls and transferred by the transition and the collisions to the fluid velocity and 
heating, and ultimately to the perturbations such as the gravitational waves. 

• The nucleation size r n of the bubbles of the transition between phases n, n — 1. 

We want to express r n in terms of other physical parameters. To proceed, we need to 
know a bit more about the process of bubble nucleation. 

The growth or not of a bubble can be seen as a competition between the expansion due 
to the release of energy from the transition and the surface tension of the bubble wal0. 
Only bubbles of an initial size larger or equal than a certain critical value can effectively 
be nucleated and grow. 

A precise description is possible in terms of a tunneling action, or free energy in case 
the temperature is not zero, for the order parameter of the transition [12]. The tunneling 
action/free energy is indeed the sum of a term relative to the bubble's wall tension minus a 
term for the energy of the bubble interior. The critical radius is computed by extremizing 
it. 

In the thin wall approximation for the tunneling, we find 

for T < r~ l (1) 

for T > r~ l , (2) 

where S n is the surface tension of the bubble wall and Ae n has been defined above. We 
have also expressed the critical radius in terms of the extremized tunneling action or 
free energy f£\ where = ^ F«> = [H]. 

The two formulas above are respectively valid for vacuum or thermal nucleatiorjf]. The 
tunneling action is used if T < r" 1 (vacuum description), otherwise the free energy must 
be employed [12]. We will show in section 13.3.11 which of the two description is more 
appropriate. 

3.2.1 Conditions on the decay rates and the time-scales of the transitions 

The phase transitions must be successful, reaching the stages of percolation, bubble collisions 
and large scale thermalization, preserving the homogeneity and isotropy of the Universe. In 
order to do so, the transition must cope with the inflationary expansion of the Universe. 

4 We neglect gravity as it will be easy to check later that r n <C Rs, where Rs ~ (GAe„)~^ is the gravitational 
radius signaling the need to consider gravity (G is newton's constant) |13j . 

5 Recall that the latter could be possible, as after the first transition completes there is also a radiation 
component of the Universe coming from the bubble walls decay. 



3^_ 
Ae n ' 

2S n (T) _ 2 
Ae n (T) (16 7T 
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This puts a series of conditions on the decay rate T n and the time scale of the transition 
[H1EH] • In particular, the necessary conditions for successful completion of a transitions during 
inflation is [H] 

f n >3H. (3) 

T n depends in general on timfl we can have fast tunneling models (as chain inflation), where 
the decay rate is large enough at the onset for the transition to complete very rapidly, and 
evolving tunneling models, where the transition occur initially with a smaller rate (the Universe 
is trapped in the old vacuum for a certain time) and finishes when the tunneling rate has 
increased enough to satisfy (J3J). 

However, to preserve homogeneity and isotropy, most of the bubbles must nucleate and 
collide in a short time, compared to the Hubble time (big bubbles are dangerous, see [H]). 
This implies 

< H" 1 . (4) 

We are now going to find more stringent bounds on f3 n considering the evolution of the 
bubbles. 

Bubbles are nucleated with radius r ~ r n . By the time of collision, a point on the surface 
of the wall of a bubble has moved by a physical distance 

r f , n -r n = - (e/»» - lj , (5) 

where v is the wall velocity. We have assumed H and v to be constant during the time of the 
transition, as the time-scale of the latter is necessarily short, as we said. Initially the bubbles 
are spherical. We assume, again because of the short evolution time-scale, that they retain that 
shape. 

In all realistic cases of successful transitions two more conditions are verified^ [SHOE] 

r f , n -r n ~^- (6) 

Pn 

r n < f . (7) 

Pn 

For fl5]) and ([6]) to be consistent at least to an acceptable value (say 5%), it must be 

|>10. (8) 
Condition ([7]), instead, gives an upper bound on (3 n . We start by writing Ae„ as 

~ ~Tt^ ~ e^PandA"^ , (9) 

where p is the total energy density dominated by the vacuum component, and we have used 
the Friedman equation. 



6 For example through the dependence of the tunneling action on different fields. 

7 In principle, we could partially relax these conditions, but those transitions would not be typical and, more 
importantly, the existing numerical studies [3H6], which we extend in section f4. 41 would not be applicable. 
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Inserting ([9]) in (JT]), we find from (J7|) 

We have considered only the case T < r" 1 as we will soon show that it is the relevant one. 

Summarizing, the bounds on the ratio between the scale of the transitions and the Hubble 
rate are 

10 <#< (£)**.. (ID 



3.3 Radiation fluid 

The collision of walls and the release of the latent heat generically produce a radiation- 
dominated fluicfi In order to have an effective Friedman-Robertson- Walker description of 
the metric, the fluid must thermalize and the time-scales for the decay of the collided walls and 
the fluid thermalization must be shortly. 

The process of bubble collision and the transfer of the energy to radiation are complex 
phenomena that can be studied numerically. We are not going to do so and we will assume 
that the thermalization of the fluid occurs rapidly. 

Note that the production of radiation cannot be neglected: although during inflation it is 
a sub-dominant fraction of the total energy density in comparison with the vacuum, it can be 
important for what concerns the perturbations. This was shown, for example, in the case of 
chain inflation with a single scalar field, where it was actually essential to provide a spectrum 
of adiabatic density (scalar) perturbations in accordance with observations [9] . 

The parameters necessary to describe the fluid are the following ones: 

• The temperature T of the radiation fluid. 

Patches of the Universe in different phases have in general different temperatures, and 
we should distinguish them with a label n. As we will see (section |4.4. ip . the differences 
between the temperatures of the phases are (and must be) negligible for inflation to be 
efficient, therefore we will often omit the suffix n. 

• Plasma scales. 

The constituents of the fluid will in general be charge carriers. Such fluid goes under the 
name of plasma. More refined definition of plasma are possible, but we will not consider 
them. 

When it is relativistic, its dynamics is regulated only by the gauge coupling(s) g and 
the temperature T. The relevant scales, such as the plasma frequency u p ~ gT or the 

8 We will show that it is relativistic (radiation) in scction l3.3.1l 

9 If it were necessary to remind it: it is well known that in all FRW cosmological models matter and radiation 
are really never in thermal equilibrium, as those space-times do not posses a time-like Killing vector. However, in 
general the Universe can be, and often is, very near thermal equilibrium and this is the meaning of thermalization 
in this context [15]. See also [9]. 
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screening distance Xd ~ (pT 1 ) l , are determined by these quantified [16]. 

Other important plasma scales and parameters will be introduced in section 14. 5[ when 
studying the gravitational emission. 



3.3.1 Conditions on the fluid temperature 

A quasi- de Sitter inflationary evolution imposes an upper bound on the fluid energy density 
Pf and therefore its temperature T. 

By using the Chaudhuri equation f!7ip and the parametrization of the Hubble constant 



H = a/ 8tt 2 c s t] e M Planck ~ 10 4 y/e M Planck , (12) 



obtained normalizing the spectrum of scalar perturbations ~ ^ — to the result of the 

«7r c s e M Planck 

5-years WMAP survey = 77 ~ 2.5 x 10~ 9 [17], we find 



T =(—Y < HTVM M Planck , (13) 



where k = k(T) counts the number of relativistic degrees of freedom at temperature T. The 
inequality follows from the fact that in the specific models there can be other components 
entering the Chaudhuri equation beside radiation, for example the kinetic energy of scalar 
field(s). 

If e is not pathologically small, T 3> GeV. At such temperature particles^] are effectively 
massless. 

Although the temperature is constrained as shown in (|13|) . we still have to check the modifi- 
cation to the tunneling dynamics due to thermal corrections. In fact, the transitions (tunneling) 
and the gravitational emission can change depending on the temperature. The vacuum tunnel- 
ing description is well-suited to capture the dynamics only provided that the temperature T is 
smaller than the inverse of the scale length of the field solution r" 1 [TJ]. From fll]), (JDJ, (JT3|>. 
we find (for k ~ 100) 

|t<( 5 #) , (4) 1 < 1 fo ^ ( "'< 20 ' < 14 > 

Noting that the critical tunneling action cannot be too large as it will suppress the decay rate 
too much, also in view of we conclude that the pure vacuum description of the tunneling 
is appropriate. 



3.4 Final comments 

The bounds on the typical scale of the transitions and the temperature (especially (ITT]) . (1TB"]) ) 
will strongly affect the gravitational emission and the sources. 

10 One generally distinguishes between different kind of collective oscillations [16], but their characteristic 
frequencies are of the same order. Also, non-abelian and abelian gauge theories are generally distinguished by 
different constant of proportionality involved in the formulas for the typical scales. We treat both cases by 
indicating the plasma scales up to the proportionality constants. 

11 We do not know the theory at those energies, it could be supersymmetric or also a GUT. 
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Let us anticipate that, in particular, the physics of the plasma will prove to be a subdominant 
source because of the constraint ( TTSj) on the temperature. The bounds (fTTj) will instead be very 
important for discussing the suppression of the waves emitted by the collisions of bubbles. 

Note also that the energy density Ae n released by a transition is much smaller than the 
total energy density, by a factor 2-g-e, see Q). This will have important consequences for the 
gravitational emission. 



4 Gravitational waves 

We organize the study of the gravitational waves in the following sections. 
4.1 Useful definitions and notation 

A gravitational wave can be seen as a ripple on the top of the background metric, as (in 
conformal time rj) 

ds 2 = a(r)) 2 {-dr) 2 + (<% + h ij )dx i dx j ] . (15) 

In particular, gravity waves are true tensor modes and are transverse, symmetric and traceless. 
The linearized Einstein equation reada^l 

+ mU l3 - V 2 hij = 167rGa 2 7rg . (16) 

where H = aH and nfj is the traceless symmetric transverse tensor part of the anisotropic 
stress tensor, obtained by suitable projection from the total energy-momentum tensor. 

Not indicating polarization indexes to avoid cluttering of formulas, the spectrum Ph(k, t) of 
gravitational waves is defined as 

(h%k,t)h(k',t)) = S^(k-k')^-P h (k,t) (17) 
and the energy density per octave 

j dp h (k,t) ^ f d „ (h*(k,t)h(k',t)) _ h3 f dQ (h*'(k,ri)h'(k',r 1 )) 



dk J (27r) 3 87rG J (2^) 3 8vrGa 2 v ' 

Here, h(k, rj) is the mode function of the graviton field h(x, rj) expanded in eigenf unctions of 
the Laplace-Beltrami differential operator with eigenvalues —k 2 . The wavenumber k is related 
to the physical momentum by k = a(t)p. 

The brackets in ( II 7p . ( II 8p refer to a quantum treatment of perturbations. We will evaluate 
the two-point functions on the Bunch-Davies vacuum, neglecting the so-called transplanckian 
issue. 



12 In cosmological perturbation theory the sources we will consider arise at second order. At that level tensor 
perturbations are gauge dependent. 
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4.2 Formal solution of the wave equation 

In order to discuss the observables we will be interested in, the field equation (1161) must be 
solved both during inflation and later on, during matter or radiation domination. 



The terms in the solution proportional to Ci,C2 are related to the vacuum fluctuations 
of the field (homogeneous equation). We choose the Bunch-Davies vacuum setting c\ = 



During matter or radiation domination, when the source is no longer active and the scale 
has reentered the horizon (k 3> Ji), the solution is 




The constants A±{k) are determined by properly matching (I19p . (|2"T|) so that both h and h! are 
continuous. The time of matching depends on the duration of inflation and on the value of krj, 
which signals when and if the scale k~ l is inside or outside the horizon. 

4.3 Sources of gravitational waves in the presence of first order tran- 
sitions during inflation 

First order phase transitions provide many possible sources for gravitational waves. In partic- 
ular we will consider 

• the anisotropy stress tensor from bubble collision 

• the velocity spectrum of the fluid 

• hydrodynamical turbulence 

• non-zero gauge fields and (hyper) magneto- hydrodynamicaf^l turbulence 

• anisotropy tensor of the radiation fluid (viscosity) 

13 We will indicate with the term "hyper" the fields associated with a U{\) gauge symmetry, although in 
general that will not be necessarily related to the hypercharge of the Standard Model. 





where H^J 2 ^ are Hankel functions and 

V 1 




H s£ e i{v T +m, C2 = . 
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In the formal expansion of cosmological perturbation theory, these sources arise at second order. 
At that order, tensor modes are also sourced by first order metric perturbations, coming from 
the expansion both of the energy momentum and of the Einstein tensors [T8]. We will not 
discuss this in the present work. 

We turn now to the detailed analysis of the sources listed above and their wave emission 
during inflation. 



4.4 From phase transitions and bubble collisions 

The collisions of bubbles source gravitational waves. Instead, the evolution of bubbles prior 
to collision does not generate gravitational waves because of the spherical symmetry of the 
bubbles [TJ. 



4.4.1 Features of the source 

We have found at the end of section l3.3.1l that the appropriate description for transitions that are 
compatible with inflation is the vacuum one, in terms of the tunneling of an order parameter 
(scalar field) [12j[T3]. There is nevertheless an amount of radiation, due to the collisions of 
bubble walls of previous transitions. 

The energy released by a transition goes partly in the acceleration of the bubble walls, and 
partly in the velocity and heating of the fluid. The features of the gravitational waves emission 
are different depending on how the energy is divided in these two channels 

We can understand how much of the energy released by the transition goes into the fluid 
by looking at the hydrodynamical equations for the fluid at the wall. In fact, the fluid can 
reach a steady-state configuration, at some time after the nucleation and at some distance from 
the wall. In that case, the transition front (related to the bubble wall) acts as a discontinuity 
surface where the energy and momentum fluxes must be conserved, leading to the equations of 
conservation in the rest frame of the walf^l 
4 4 

3 KT n V (f),n = 3 KT n-l V (f),n-1 (22) 
-kT^ ^(/) )n + g K ^n — e n = v ff),n-l + g^n-l ~~ e n~l ■ (23) 

Here W(/), n (n-i) is the component of the four- velocity of the fluid in phase n{n — 1) locally 
orthogonal to the wall. 

We try to find what can be the values and solutions for the velocity that comply with 
inflation. Considering that inflation allows only small velocity and temperature perturbations, 
$22]), (|23D lead to 

Ae^ < l \8(Klt_ x )\ 

nT± ~ 3 kT± ' 1 ] 



where 8(kT^_j) = kT^_ x — kT^. 
But from (jS]), (17T|), we find 



14 Hcre, for the moment, we have assumed that the interactions between the fluid and the wall are negligible 
during the bubble expansion prior to collisions, but it is possible to extend this to the more general case (see 
for example [T§]). 
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whikH 

l^(^ 4 -i)l = M w j^Aoti ^ g (26) 

«^ pf ptot Mpi anckv /cl?' 

By looking at fllip . (IT2"|) . we see that (I25p . (125]) are in contradiction with ([Ml); which must 
be satisfied if there exist steady state hydrodynamical solutions at the wall complying with 
inflation. Therefore there is no such solution. 

The meaning of this result is that, for successful transitions compatible with inflation, the 
vacuum energy released in the transition goes mostly into the acceleration of the walls, which 
expand at a velocity rapidly approaching the one of light until collision. We can therefore assume 
v ~ 1 for the wall velocity. The bubble dynamics is well described by the approximation of 
bubbles in vacuum. 

The gravitational waves are therefore mainly sourced by the stress tensor of the scalar field 
describing the bubbles' configuration in vacuum, with a sub-dominant contribution from the 
fluid, suppressed by powers of its small velocity. 



4.4.2 Calculation of the wave spectrum 

We turn now to the detailed computation of the spectrum of gravitational waves. They are 
sourced because of the breaking of the spherical symmetry of bubbles at collisions. The latter 
are complicated events that occur out of equilibrium. It is therefore quite difficult to give an 
analytic description of their gravitational emission: in principle we could expect a whole range 
of different scales to appear in complicated inter-correlation. 

The only way to deal with these problems is via numerical simulations [3HB] . On top of those, 
useful analytical formulas can be derived [S1E]- We will obtain the spectrum of gravitational 
waves by extending the results of the numerical simulations to the case of transitions occurring 
during inflation. 

The best simulations available today [3H6] consider short-lasting sources and static back- 
ground (neglected expansion of the Universe, no inflation). For nearly vacuum collisions and 
wall velocity v ~ 1 in their formulas, they show that only two quantities are important in 
determining the spectrum of the gravitational waves: the overall energy density pw released 
by the transition and carried by the walls that collide, and the time-scale of the phase 
transition, which sets the peak frequency of emission 

References |5,6J also show that the energy density per octave radiated by colliding bubbles 
goes like k 3 for small frequencies, while for large ones it decays as k" 1 . In particular, there 
is a single peak when k ~ f3 n and the maximum (normalized to the total energy density) i^ 

0.013 pw 



15 We estimate the enthalpy density variation compatible with inflation using the value on supcrhorizon scales 
as it is larger than that at sub-horizon scales. The result for density perturbations that we use in the formula 
can be rigorously obtained in the flat gauge. We consider models of inflation, where isocurvature perturbations 
are negligible. 

16 Hcre. we specifically use the results in the more recent reference [5]. In adapting from their conven- 
tions to ours, one has to compare the formula we obtain using (TTO)) in |[T5]). for large k, with Weinberg's 
formula from |20j . used in [5]. Note that we define the Fourier transform of the energy-momentum tensor as 
Tfj, v (k,k) = J (2tt)^ elk x Tij,v(x, t), whereas [5] uses Weinberg's normalization for the Fourier transform. The 
physical equivalence of Weinberg's formula and the one we obtain from (fT5)) . (flT))) for large k has been discussed 
ini. 
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We stress that the transitions considered by [3H6] occur at the end or after inflation and the 
energy carried by the walls in their cases is nearly all of the total energy, so that 8nG pw = 3iJ 2 . 

We will now derive the radiated energy and spectrum of waves in our scenario. There are 
three evident differences with the setup of [3H6]: 

• the Universe in our case undergoes inflation. However, the static approximation at the 
time of emission of the waves can be acceptable as long as their modes are k > aH. This is 
compatible with our scenario as the typical scales/frequencies /3 n of the phase transitions 
are larger than the horizon (see section 13.2.11) . 

• in our case, the total energy released by the walls is Ae n given by 0, which is much 
smaller than the total vacuum energy driving inflation. 

• the evolution of the waves during inflation is peculiar as modes can exit the horizon. 

Using this information, we can now adapt the results of the simulations and derive the energy 
density per octave radiated by colliding bubbles of the n-th transition during inflation. We 
find that initially, at the time t n when the sources stops being active shortly after the collisions 
(time of emission), 



Ptot dp 



t=U. 




H n <p n < p n 

Pn > 



where p tot is the total energy density. We have written the result in terms of the physical 
momentum p n = ka' 1 , a n = a(t n ). 

When k 3> "H, there is a simple relation between the spectrum (fiTj) and the energy density 

ak Sna^G 

Using it, we obtain the spectrum of gravitational waves at the time of emission t n 

(29) 



t — £71 




We need now to evolve the spectrum after the moment of emission of the waves. As this 
evolution is particular during inflation, further differences with respect to the results in [3H6] 
will arise. 

At the time of emission (roughly at the completion of the transition and collisions), the 
modes are within the horizon. From f[T9~j) . we can see that then the wave evolves as ~ a -1 
until the mode exits from the horizon at time £ ex . After that, the wave remains constant until 
reentering horizon. 

We need to compute the proper redshift of the modes at t > t n . For modes still within the 
horizon at time t, we find = — . For modes that have already exit the horizon at some 

' a(t) p n J 

time t cx < t, instead, the redshift is only — = — —. Overall, the redshift factor at time t for a 
momentum p can be written as 

n p n j j- superhorizon modes at time t {pit) < H) 

— Xn j Xn = \ n , , , , . ,//,\ rT i \^) 



a(t) pn' }-£- subhorizon modes at time t (p(t) > H) . 
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where we have used H ex ~ H ~ H n , since the Hubble parameter does not evolve much during 
inflation. 

Therefore, at time t > t n (still during inflation), the spectrum generated by the collisions 
of bubble of the n-th transition has evolved to 

p > -^{wj ex '\(^y P > /A , ■ (31) 

where 

fPn = 77\~ Pn > /ff„ = —7JT- Hn ■ (32) 

a(t) a[t) 

If we want to have a chance to detect these primordial gravitational waves, we must consider 
super horizon modes, as sub horizon ones are more suppressed by the redshift. Their total 
spectrum at time t is obtained by summing over all the phase transitions and is given by, using 
( 130|) . ( )3~Tj) for superhorizon modes, 




where 



a(t n )H < k < f3 n 



-^Planck \ 2 ( H 



(33) 



/3 n = a n p n , k = a nPn >a(t)H, A n = 0.07\^^j (-1 e 2 (34) 

This represents the main novel result of this section. 

We will deal with the later evolution of the modes, after reentering the horizon, in section 

E2 



4.4.3 Comments on the final result 

Comparing the spectrum (133]) with the one generated in vacuum during inflation (see formula 
(|82|) in appendix IA.3j) . we see that the main differences follows: 

• the spectrum (1331) is not scale invariant 

• due to ffTll . ffT2l) . A n is an additional suppression factor, with respect to the waves 
generated by vacuum oscillations 

• the accumulation due to a sufficiently large number of phase transitions N could cope 
with the suppression and make the waves detectable. 

In discussing the backreaction of the waves (and possibly the breaking down of the perturbative 
approach), we need to consider the possible number iV of transitions. 

As we see, for a few phase transitions, the waves sourced by the collisions of bubbles are 
even more suppressed than those produced by vacuum quantum oscillations, and therefore their 
backreaction is negligible (perturbation theory is accurate). 
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For a large number of transitions, we have in principle a stronger backreaction. It is easily 
verified, though, that the number of transitions needed to affect the background evolution is 
unlikely to occur (for the range of parameters given in section [3] and e ~ 10~ 2 , it would have 
to be at least of the order of N ~ 10 12 ). 

We might also wonder if the depletion of the energy density due to gravitational radiation 
could affect the scalar density perturbations. However, from ( 133|) and the above comments, it is 
evident that this depletion of energy for a single transition is very small (smaller than the linear 
order in perturbation theory). Only a large number of transitions could have an effect (in that 
case, invalidating the perturbation expansion). However, even in that case this would happen 
after many transitions have occurred, that is quite late during inflation. It would therefore 
affect only the late density perturbations that are not interesting for detection at the CMBR 
(see also section 15. ip . 

The spectrum ( 13 3 p is different also from the one generated by the collisions of bubbles of 
first order phase transitions that do not occur during inflation. Indeed, beside the fact that 
in this second case the fluid velocity spectrum is generally an important source and that the 
maximum amplitude of the waves from a single transition is less suppressed, the most notable 
difference is that in our case the presence of modes that were superhorizon at inflation gives 
the spectrum a different dependence on the frequency (compare with [3H7]). If detected, this 
would help in distinguishing the two cases. 

Another comment concerns the sum we have performed in fl33l . In fact, there is a non-zero 
probability of nucleation of bubbles of phase, say, n — 2 when the bubbles of phase n — 1 have 
not yet collided. Would this make the simple sum rule we use inadequate? The answer to 
this question is that it would not, as long as the moments of collisions of consecutive phase 
transitions are spaced out. Indeed, the evolution of the transition before collision (expansion of 
the bubbles) does not generate gravitational waves, if the bubbles are spherically symmetric [7]. 

4.5 From turbulence, (hyper) magnetic fields, viscosity 

The presence of a radiation dominated fluid is a consequence of the processes of bubble collisions 
and walls decay after each phase transition. Although sub-dominant, this component cannot 
be neglected when discussing perturbations, as it is can be even essential for having the right 
spectrum of adiabatic density perturbations, for example in chain inflation [9] . 

The physics related to the plasma or hydro- dynamics could generate gravitational waves 
during inflation, depending on the gauge coupling and the temperature of the fluic0- This is 
what we want to investigate now. 

We will make the assumption of weakly coupled plasmas, because their stress anisotropy 
tensor is larger (large viscosities) than that of strongly coupled plasmas, and therefore a po- 
tentially stronger source of gravitational waves. Nevertheless, the anisotropy tensor cannot be 
too large, otherwise it would spoil the homogeneous and isotropic description of the Universe 
at large scales. It will appear that this indeed does not happen when the parameters satisfy 
the bounds in section [3j 

Note also that the validity of the hydrodynamical or plasma description depends on the 
relative magnitude of the Hubble scale and the microscopical ones, such as the mean-free-path 
or the screening distance. This tells us that only for some ranges of values of the couplings 

17 For an incomplete bibliography on the processes themselves and the associated gravitational production 
(not during inflation), see [2T| . 
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the description is self-consistent. In particular, the couplings cannot be too small, although a 
precise bound depends on the model-dependent numerical factors entering the formulas of the 
plasma scales. 

We will make use of the transport coefficients characterizing the transport of energy, mo- 
mentum and charge across the fluid [16]. They have been generally calculated within the most 
well-known theories (especially QED, QCD), but usually the results can be extended with 
only minor modifications to different higher energy theories, thanks indeed to the plasma or 
hydro dynamical approximation. We assume that this is our casd^l. 

We are now going to analyze the various sources related to hydro- and plasma dynamics. It 
will turn out that they are very much suppressed, contrarily to what generally happens when 
first order transitions occur not during inflation. The ratio J- will be particularly important in 
the following considerations. It is strongly constrained by the requirements of efficient inflation 
and self-consistency of the theory (see section [3]) and this will ultimately be the reason why 
these sources do not produce a sizable spectrum of gravitational waves in our scenario. 



4.5.1 Hydro dynamical turbulence 

Turbulence is a strong source of gravitational waves. Let us study its occurrence in our scenario. 
In our case, the generating mechanism would be the collisions of bubbles, stirring the fluid. The 
typical length scale of injection is approximately the bubble size at collisions, ~ f3~ 1 , and this, 
together with the kinematic viscosity [16] 



^g-Hogig- 1 )-^- 1 , (35) 



yields the Reynold number 



Re = ~ vf^W logGT X ) , (36) 

V Pn 

where vf" ] is the characteristic velocity of the fluid flow at the injection scale. 

Large Reynold numbers signal the onset of turbulence. We therefore see from (1361) that the 
relevant condition is 

T 1^1 turbulence . . 

£ = (<! no turbulence. (37) 

Looking at (TT2} [TBI [TT]) , we find that in our scenario 

^<4. (38) 

Pn 

Therefore, turbulence does not occur and a sizable emission of gravitational waves is not possi- 
ble. The reasons for this are the smallness of the scales of injection of energy and the rapidity 
of the bubble evolution, which are consequences of the requirements of small backreaction on 
inflation and self-consistency of the theory (section [3]). 



18 Our results will therefore apply to theories and models for which this is possible. In the following, our 
formulas for the parameters will be written up to proportionality constant (indicated with a index o), which 
depend on the peculiar details of the plasma under consideration (such as the number of light versus heavy 
species or the rank of the gauge symmetry group). These will not be influential for the results. 
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4.5.2 Plasma physics and gauge fields 

The radiation fluid is likely charged (therefore it is generically a plasma) and local charge 
asymmetry and currents can be generated^!. We investigate now whether non-zero long-range 
gauge fields generated by these currents could represent a sizable source of gravitational waves 
during inflation. 

The only long-range fields that can survive in a plasma for enough time are magnetic 
fields associated with a £7(1) gauge symmetry (hypermagnetic fields), see section 5.6 in [25] . 
The precise form of the generated fields depends on the details of the bubble growth, the 
couplings, the tunneling events, the non-equilibrium dynamics. Nevertheless, it is possible to 
draw sufficient conclusions from general considerations. 

The fields are generated by the currents at the scale L p ~ [3~ l of the bubbles^ (21] and 
then possibly enhanced by magnetic turbulence, which is signaled by a large magnetic Reynold 
number. This is defined as 

Re, = , (39) 

where L p is the typical scale at which magnetic fields are generated and a c = °"o g 2 lo ^ g -i) is the 
(hyper)electrical conductivity at high temperature. 
In our case, 

Re »~ V f n) -R J ( IV ( 40 ) 

which, even for reasonably small couplings (g ~ 0.1) is much less than 100, given the bounds 
in section [3J Therefore there is no hypermagnetic turbulence and it is unlikely that the small 
scale hypermagnetic fields appreciably source gravitational waveJ^I. 

4.5.3 Stress tensor from fluid viscosity 

The form of the anisotropic stress tensor in general relativity for a fluid with velocity u is [20] 



< isc) = -? (d jU P + d iU f - 2 -V ■ ti<%) - CV • tf<%. (41) 



C is the bulk viscosity. For a relativistic fluid, £ is vanishing, as it is given by 

i dp f 



C = CoTV 



3 dP f 

where p/, P/,n/ are respectively the energy, pressure and number density. 



(42) 



19 For example, it is enough to have different mean free paths (different couplings) for the various species to 
generate a local charge asymmetry, in presence of particle-antiparticle number asymmetry. Then, the collisions 
among bubbles impart a vorticity to the charged fluid, with the creation of microscopic currents. The num- 
ber asymmetry can occur because bubble collisions and first order transitions are likely to provide the right 
conditions 22 . Nevertheless, a precise statement can be made only with a detailed model pH] , 

20 Due to the different orientation of field lines at the bubble scale, it is also possible to assume that the field 
spectrum is stochastic, in this way preserving the global isotropy of spacetime. 

21 In principle, there is still a residual but unlikely possibility that in specific models even the small scale 
hypermagnetic fields can be strong enough to generate a sufficiently sizable amount of gravitational emission. 
We will not discuss this case. 
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The part of 7r|-J lsc ' 1 proportional to the shear viscosity q could instead be important, as 

T 3 

q = ^O^j 7 7T (43) 

for a weakly coupled fluid. Nevertheless, it turns out that this source is negligible as well. 

In fact, the energy per octave radiated by an anisotropic stress tensor 11^, normalized by 
the total energy density, at emission time 77 = on a scale k > H is 



_ 1 dp h {k) 



Ptot dk 



16 f 

Gk'al / dQ{U. ij (k, fc) n *#, k)) , (44) 

Ptot J 



where a* = a(r]*), while Uij(k, k) = J e lkv lk ' s Ilij(x,r]) and / dQ is the angular integral. 
In the case of (jHJ), considering only the shear viscosity part, we have 

fi (visc) ^ 1^! V 2 (45) 
Ptot 

where V 2 = alk 5 \ 5 Vg A and Vg A (AA;, k) is the dimensionless square mean-field gradient of the 
velocity (integrated over the angles and Fourier transformed in the time dependence) obtained 
by smoothing over the comoving length scale A and in which we have also absorbed the de- 
pendence on kX from the velocity spectrum. The scales suitable for hydrodynamical treatment 
must be larger than the coherence scale of the velocity, therefore in our scenario A ~ /3" 1 . 

On the other hand, at the moment of emission, the energy per octave radiated by the 
collision of bubbles of one phase transition is (1271) . here written for comoving wavenumbers, 




fl M) = ^ £L £ 2 [Pn\ (46) 



We have considered wavenumbers appropriate for the comparison with the scales of the hydro- 
dynamical description, which means k < j3 n . 

For scales k < A^ 1 , where A/ = Aog -4 log((? _1 ) _1 T _1 is the mean-free-path, we therefore 
find 

Q(visc) 

T^Hy <v 2 9 J k=xjl <l, (47) 

where we have used (TTT1 [T2l [T3|) and the final estimate takes into account that v\ A is a small 
perturbation during inflation. 

The energy radiated in gravitational waves from viscosity is therefore smaller than the one 
coming from bubble collisions. 



4.5.4 Thermal fluctuations 

There are also other possibilities for generating (hyper)magnetic fields in a plasma, beside local 
currents from bubble dynamics: for example quantum and thermal fluctuations. In both cases 
the generated fields can source gravitational waves. 

A complete analysis of these sources would go beyond the scope of this paper, as they are 
not strictly present only when first order transitions occur and the results would depend on 
specific details of the models. 
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In fact, quantum fluctuations can generate gauge fields extended at appreciable scales only 
when their couplings to the metric are not conformal symmetric (see for example [21] and 
reference therein). On the other hand, the excitation of gauge fields via thermal fluctuations has 
been studied in flat Minkowski space via thermal field theory or the Boltzmann equation [25] . 
and the extension of those results to a rapidly expanding Universe is not straightforward. 



5 Detection 

Our analysis indicates that the dominant contribution to the spectrum of gravitational waves 
from first order phase transitions during inflation, is sourced by collisions of bubbles and has 
the spectrum (133|) . 

We will now discuss the possibility of detection of this spectrum both in the CMBR 
anisotropies and by direct measurement at interferometers. Although detailed, our analysis 
will not be precise down to numerical factors of order one, which depend on the particular 
high-energy models. 



5.1 CMBR 

Primordial gravitational waves affect both the temperature anisotropies and the polarization 
of the CMBR. Experiments such as CMBPol and Planck will investigate these observables in 
the next future. Here, we will concentrate on the temperature anisotropies, considering the 
so-called tensor (T)-to-scalar (S) ratio 

r = Zf' (48) 

where Cf, Cf come from the decomposition of the spectrum of temperature anisotropies in 
two directions li, 1 2 with l x • 1 2 = cos(6 l ), 

^(li)f?(l 2 )) = ^ f> + l)(Cf + C$)P e (co S 9) . (49) 

' 1=2 



The precise computation of r would require a numerical evaluation of (]49p using the spec- 
trum However, for our purposes it suffices an approximate analytical computation. We 
follow [26] , assuming t > 10. 

Using k = ^jr, do = 1, we obtain 

i(£+l)Cf -^Pr^f 2 ) (50) 



and therefore, for the spectrum (1331) and £(£ + 1) Cf 



4 H 2 

25 4nc s eM^ nck - 



N 



22 Recall that the multipoles we can consider arc such that I < 200. 
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where we have considered only wavenumbers in the range a n H < k < j3 n since the contribution 
from larger wavenumbers is very suppressed (see (1331) ). 

Equation f[5"Tj) is very different from the result for tensor modes generated by vacuum fluc- 
tuations in single-field slow- roll/chaotic models, which would be r = 8e: indeed, (l5"Tj) shows i) 
a dependence on the multipole t due to the lack of scale invariance in (133|) . ii) an additional 
suppression factor A n (l34l) . and in) a possible accumulation due to many phase transitions. 

We would like to understand how many transitions are necessary to have a detectable 
amount of tensor modes (we require r ~ 0.07). To proceed, we will make a simplifying assump- 
tion: a negligible dependence on n for the ratio Jr-. We will also take e ~ 10~ 2 as a typical 
value, for illustrative purposes. 

We also limit our estimate to the multipoles corresponding to the frequencies p = (3 n , 
which is the typical one of the collisions, and p ~ H n , which is the lower limit of validity for 
the numerical simulations we have extended. We indicate these multipoles in the formulas as 
^* = {^p n ^H„}- F° r these specific frequencies, also because of the simplifying assumption we 

have made above, the ratio is independent of n. 
We define 

\A n l = A-f ~ 0.07ATL0V 1 ( tL ) £ 2 ^, 

71=1 



0r, 



where we have used ( IT2|) . (134]) and = J^-. 
Taking into account (fTTI) . A lies in the range 



1.1 x lQ-^Ne ( % ) <A<0.07Ns (53) 



7T 

and therefore the scalar to tensor ratio (15TD falls within the interval 



-5 W J / ^V ^ / „ / n n„ at .2 h 



1.5c s 10" 5 Ne z [ ^ -f < r < 0.9c s N e 2 . (54) 

If k = f3, we see that TV ~ 7.7 x 10 2 transitions with the minimal allowed transition rate 
(jj ~ 10) would be sufficient to have r ~ 0.07. But the number of transitions necessary for 
detection rapidly increases for faster rates, up to TV ~ 3 x 10 8 for the fastest rate in (fTTI) . 



If, instead, it is k — a n H n , we find that iV > 77 transitions with slow rate (-^ ~ 10) would 
be detectable via the CMBR, but for the fastest transitions we need at least ~ 6 x 10 6 of 
them. 

It appears therefore that a not too large number of transitions occurring at a slower rate 
could be detectable, especially for lower multipoles. However, as the rate of the transitions 
increases, the spectrum rapidly decreases and becomes undetectable. 

Note that in many of the studied models of transitions outside inflation it is found that the 
ratio jj is larger than ~ 10. However, the situation for first order transitions in models and 
theories valid at the scales of inflation could be different. Because of our limited knowledge of 
the relevant theories at very high energy, such as the string landscape, we cannot discard those 
values of It is therefore important to analyze this point in future research within concrete 
models of interest. 
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The suppression of the emission from fast transitions also occur for transitions not during 
inflation, but in the inflationary scenario the phenomenon is much more pronounced as the 

suppression factor goes like (jj^j e 2 for each transition, see fl33|) . 
5.2 Direct detection 

We will investigate the direct detectability of the gravitational signal by interferometers con- 
sidering the cases of LIGO, LISA and DECIGdFl To do so, we need to evolve the spectrum 
after inflation ends, from the moment when the mode reenters the horizon until now. 

Since, differently from the previous sections, we are now evolving the waves also after infla- 
tion, we will indicate as H in Q the Hubble rate during inflation, to avoid any possible confusion. 

From (12 ip . we see that after reentering the horizon, h evolves as ~ a -1 . We define the 
transfer function Tip) = where a p is the scale factor at the time of reentering for the 
physical momentum p measured today. As discussed in \27\ . the sensitivities of LIGO, LISA, 
and DECIGO peak around frequencies which had to be within the horizon well before matter- 
radiation equality and nucleosynthesis. Therefore, assuming adiabatic expansion after the end 
of inflation, [27] 

rr-i \ 9n /0.63Hz\ /1(D0\ 5 

Tip) = 2.1 x 1(T 20 (55) 



V ) 

where k p is the number of effective relativistic degrees of freedom contributing to the energy 
density at the moment of the reentry of the scale p~ l . 

We discuss the possible detection of the primordial waves in terms of the strain amplitud^ 



^ = T(ph/^M. (56) 



By using (|5"5"|) and we find 



V / \ «p / „ V^Mpianck V V V P 

where we have considered only wave-numbers a n H m f{ < k < j3 n , since the contribution from 
larger ones is very suppressed, and we sum over the different transitions. 
Here, 

a n (3 n n v 14 /100\MGeV 



/A.,0 = =^ = ^-8.0 x 10-" — (58) 

a \ ^end / 1 end 

M n is the number of e-foldings from the moment of collision of bubbles of the n — > n — 1 
transition until the end of inflation, and T en d, ftend are evaluated at the end of inflation. In the 
following we will w 10 2 for illustrative reasons. These are also the typical 

values in GUTs and minimal supersymmetric models. 



23 We will actually consider Ultimate DECIGO. 

24 0btained re-writing formula (3.2) in [35] with our conventions in (jTTJ) , (p~5|) . See also section 2.2 in 
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To evaluate (jSTl l58|) . we need to compute the temperature T en d, which is given by 



JV 



T cnA = Y,T n e- N - + T f . (59) 



n=l 

Here T n ~ T is the temperature of a single phase transition, see ffT3]) . and 7/ comes from a 
possible final decay of the inflaton (reheating), which could be also one last phase transition 
with large backreaction. 

We can perform the sum in (1591) by considering that the last transition occurred A/i ast e- 
foldings before the end of inflation and that the earlier transitions and collisions times were 
spaced out by intervals of AAf n ~ AAf e-foldingsEl. Then we define 

= \ _ e e _^ = F(N, AAf, Mast) ■ (60) 

n=l 

Assuming finally, for simplicity, that Tj < ^ n ^ 1 n e _A/ ' 1 , we obtain 

Tend « T F(N, AAf , A/iast) • (61) 

Let us recall now that (1331) is valid only for p n > H n . These high frequencies could fall 
within the range of sensitivity of the interferometers only after sufficient redshift. The sum 
over transitions in (]57|) will therefore start from the transition n m i n for which at least the 
smallest frequency that we can consider at the time of production t nmin has had the necessary 
redshift until now. That frequency is p Hn ™in (t„ min ) = H ninin ~ H in Q, which is redshifted today 
to ^ 

p »^ M = ^p--...(f„ ) = 1.2 x 10^ F(jv ^'^ t) H Z (62) 

Comparison with experimental setups 

We are now ready to compare our results with the sensitivities of LIGO, LISA and Ultimate 
DECIGO. The latter peak around certain frequencies as follows [28|l29]: 

(63) 
(64) 
(65) 

The results we list in table [T]have been obtained by asking for the strain amplitude fl57|) to be 
within the sensitivity of the detectors when evaluated at the respective peak frequencies. The 
strain is computed using flBTt I58| IBT] I6"i2"|) . We have chosen ~ 0(1), e ~ 10^" 2 as indicative 
values and assumed for -|r- a negligible dependence on n. The range of values considered for 
A = J2 n A n is given by ([53]). 

The detectability of the emitted waves depends on the number of transitions, their timescales 
and the number of e-foldings that must have been occurred to sufficiently redshift their fre- 
quencies after the waves were emitted. These quantities enter the strain through the function 



LIGO : 


h, 


~ 10~ 23 Hz"5 


at 


/~ 


100 Hz 


LISA : 


h. 


~ 4 x 10" 21 Hz"^ 


at 


/~ 


10" 3 Hz 


UDECIGO : 


h, 


~ 10" 27 Hz"5 


at 


/~ 


0.1 Hz. 
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f (H7) 
J 


Detectability 






No if 


"Yp^i for all no^iTilp — if 

-L \yiD L\JL CLL1 pUDOlUlC rj 1-L 


t tpo 


1UU 


F(N-n min , *f,0) < 1.6 x 10 6 


F(iV-n m in, > 4.4 x 10 8 


LISA 


10~ 3 


F(N-n min , M:,0) < 19.7 


F(N-n min ,^f,0) > 5.5 x 10 3 


UDECIGO 


1 


F(iV-n m in, < 4.9x 10- 3 


F(N-n min , *f,0) > 1.4 



Table 1: Bounds on N,n min , AAf via the function F for the direct detection of primordial 
gravitational waves from phase transitions during inflation. 



F(N-n m i a , 4^,0) (see (1601 ). Table [T] shows the detectability bounds on this function for given 
N,n min ,AM. 

To understand better what these results mean for the physical parameters, we specialize to 
the case of chain inflation [9] as an illustrative example. This one is a fast tunneling model, 
where inflation stops shortly after one last phase transition (so that A/last ~ 0) and where 
N > 1 and AAf ~ Sip. From HQ]), in this case 



T?( AT A ^ ) 2 lh fOT N ~ nmin » /i' A » 1 ^ 

F(N-n mm , ,0)~< mfl R mfl !, nfl . (66) 

1 2 J \N-n mhl for /-»A-r imin , » 1 

Using f )66|) . f TTTT) and table [H we see that, for instance, in the case N — n m j n ^> jr— (many 
transitions), the gravitational waves produced by bubble collisions in chain inflation will be 



• detectable at Ultimate DECIGO for all allowed values of 

• detectable at LISA for jf- < 11. 

In particular, if we take the most favorable casd^l. j^— ~ 10, the detected waves would have 
been emitted by phase transitions occurring at least 20 e-foldings before the end of inflation 
for Ultimate DECIGO, or 27 for LISA. Unfortunately, the gravitational waves would be unde- 
tectable at LIGO. 

If detected, these waves could be distinguished from those of transitions occurring outside 
inflation, thanks to the different frequency dependence of the spectrum of modes that were 
superhorizon at inflatiorJ^I. 



25 Here, we have made the simplifying assumption that the intervals of c-foldings are approximately the same 
for all transition. The crudcness of such an approximation can be determined only by building detailed models, 
but this goes beyond the scope of this paper. 

26 Recall the comments at the end of section RTTl 

27 Note however that the strain amplitude for these waves decays more rapidly with the frequency: as p~ 2 for 
frequencies smaller than the redshifted scale of the transition and as p~ A for larger frequencies. The signal from 
transitions not during inflation goes instead like p _1 and p~ 3 respectively in the two ranges of frequencies. 
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5.3 Nucleosynthesis bound 

There is an important constraint on gravitational emission: the waves that reentered the horizon 
before nucleosynthesis must not interfere with it, and therefore satisfy [28] 

1Hz 

h p < 1.99 x 1(T 21 — . (67) 

We study this constraint just in the case of chain inflation, as an example. From formulas (15 7j) . 
(}6"6"]) . we see that in that case f lBTl) implies 

P 



Hi 



mil 



>0.21 Vp, (68) 



which is certainly satisfied by the range (I lip of allowed values for jf— ■ We conclude that the 
bound from nucleosynthesis does not rule out chain inflation. 



6 Conclusions 

In this work we have studied the production, features and detectability of gravitational waves 
in models of the early Universe where first order phase transitions occur during inflation. 

We have described these scenarios via some physical parameters, whose values have been 
constrained and bounded by an analysis of the self-consistency of the theory (in particular 
efficient inflation taking place, homogeneity and isotropy at large scales). 

The emission and features of gravitational waves are strongly affected by these bounds and 
by the specific dynamics during inflation (such as the exit from horizon of the modes). The 
resulting spectrum is different from the one due to vacuum oscillations or first order phase 
transitions occurring not during inflation. 

The first important feature is that the waves from a single transition during inflation are 
very much suppressed but the accumulation due to many transitions could make them sizable. 
Second, turbulence and (hyper) magnetic fields are a negligible source of waves, contrary to 
what generally happens when the transitions occur not during inflation. The collisions of the 
bubbles at the end of the transitions represent the prominent source of waves, yielding a non- 
scale-invariant spectrum with different frequency dependence for modes that exit horizon during 
inflation, compared to the spectrum sourced by transitions not occurring during inflation. 

We have also studied the experimental detectability of the waves. The main points emerging 
from this part of the analysis are that: 

• a not too large number of slow (but still successful) transitions occurring during the 
inflationary era could leave observable marks in the CMBR anisotropics (large tensor-to- 
scalar ratio for reasonable values of the parameters). As the fastness of the transitions 
increases, the signal rapidly weakens, requiring accumulation from a large number of them 
to be measurable, 

• direct detection via interferometers could be possible at LISA and Utimate DECIGO for 
modes that were superhorizon at inflation and for a large number of transitions. However, 
at LISA the detection could occur only for the most optimistic scenario (slow transitions). 

We could distinguish transitions during inflation from those outside inflation thanks to 
the different frequency dependence of the spectra, 
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the nucleosynthesis bound is easily satisfied by the models, in particular chain inflation 
is not ruled out. 



A Appendices 

A.l Friedman and Chaudhuri equations 

With a Friedman-Robertson- Walker ansatz 

ds 2 = -dt 2 + a(t) 2 dx 2 (69) 
for the background metric, the Friedman and Chaudhuri equations read 

H = -4nGY,(pe + Pi), (71) 

t 

where p£, Pi are respectively the energy and pressure density of the component i of the Universe 
(in our case scalar fields, radiation). The sum is over all components. Note that during inflation 

Ptot ~ /-Vacuum- 

The conformal time is related to the cosmic one by dr\ = 



A. 2 Timescale of transitions 

We call p n (t) the vacuum persistence probability, that is the probability for a point in the 
Universe to remain in the n-th vacuum at time t. By neglecting the possibility to tunnel 
directly to distant phases, it obeys the most general equation 

Pn = ~r n p n + r n+ i p n+ i (72) 

The conclusion of the phase transition is indicated by several markers [HUH]: i) the time 
t Ctn , when the patch of the Universe occupied by the old phase starts to contract, ii) the time 
tp tU , when percolation occurs, Hi) the time t s>n , when the probability p n (t) for a point to remain 
in phase n at time t has dropped below a suitable small number. 

This last requirement is actually not a faithful signal of the completion of the phase transi- 
tion (recall for example the issues in Old Inflation), but it yields an indication that is generally 
in accordance with the other two more significant conditions, when they occur, beside being 
physically reasonable. For a rapid transition it can be shown that t S)U > t c ^ n ,t p ^ n . 

The time-scale of the transition can be defined as [TJ] 

Pn 1 ~ ts,n - U,n (73) 

where t s>n ,ti in are such that@ 

Pn(ts,n) = e~ M < 1 Pn+l{ts,n) = e~ M ~ q < 1 (74) 

Pn(ti,n) = e~ m ~ 1 Pn+l(ti,n) = e~ M+q ' « 1 (75) 

28 Note that also the earlier and successive phases are sub-dominant at £s,nj but in order to estimate /3" 1 we 
need to consider only the phases n and n + 1. 
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for suitable M, M — q, m + q' 3> 1, m < 1. 

If we have T n > H from the onset, and therefore the phase transitions are occurring very 
rapidly, we can expand 

Pn(ts,n) ~ Pn(U,n) + Pn(U,n)(t s ,n ~ U,n) (76) 

and from (1721). we obtain 



P„.(^n)-p„.(^„.) = e- M -e- m (77) 
(-f n (ts,n ~ U, n ) + l) e~ m = (-f n+1 e 9 ' (t., n - t i>n ) + l) e" M . (78) 



Since e~ M < 1 



^s,n ^i,n ~ r ?l — Pn • (79) 



For tunneling rates depending on time, expanding the tunneling action/free energy around 
t Stn as S^\t) ~ (t Sjn ) — f3 n (t — i s ,n), one can also find for the decay rate per unit time and 
volume IT] 



(n) 



50) 



As we see from (f?9| [80]) . /3 n is therefore directly related to the fundamental physics governed 
by the tunneling action. 

Finally, the time-scale of the transition in terms of conformal time is 

p- 1 = Vs,n-Vi,n = a(t n )- 1 (3- 1 . (81) 
A. 3 Quantum vacuum fluctuations 

As a useful reference and comparison, we report here the spectrum of the gravitational waves 
generated via fluctuations in vacuum of the gravity field during inflation. The wave equation 
is the homogeneous version of ffTBI) . with irjj = 0. Its solution leads to the spectrum for 
super horizon modes 

p — IG£J (82) 

using the Bunch-Davies vacuum. H is the Hubble parameter during inflation. 
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